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Abstract 

We calculate the gradient of the radiation field generated by a polarization current with a 
superluminally rotating distribution pattern and show that the absolute value of this gradient 
increases as Ft 1 1 2 with distance R, within the sharply focused subbeams that constitute the overall 
radiation beam from such a source. In addition to supporting the earlier finding that the azimuthal 
and polar widths of these subbeams become narrower (as -R -3 and R , respectively), with distance 
from the source, this result implies that the boundary contribution to the solution of the wave 
equation governing the radiation field does not always vanish in the limit where the boundary 
tends to infinity (as is commonly assumed in textbooks and the published literature). There 
is a fundamental difference between the classical expression for the retarded potential and the 
corresponding retarded solution of the wave equation that governs the electromagnetic field: while 
the boundary contribution to the retarded solution for the potential can always be rendered equal 
to zero by means of a gauge transformation that preserves the Lorenz condition, the boundary 
contribution to the retarded solution of the wave equation for the field may be neglected only if 
it diminishes with distance faster than the contribution of the source density in the far zone. In 
the case of a rotating superluminal source, however, the boundary term in the retarded solution 
for the field is by a factor of the order of R 1 / 2 larger than the source term of this solution, in the 
limit where the boundary tends to infinity. This result is consistent with the prediction of the 
retarded potential that the radiation field generated by a rotating superluminal source decays as 
R~ l l 2 , instead of -R -1 , and explains why an argument based on the solution of the wave equation 
governing the field in which the boundary term is neglected (such as that presented by J. H. 
Hannay) misses the nonspherical decay of the field. Given that the distribution of the radiation 
field of an accelerated superluminal source in the far zone is not known a priori, to be prescribed as 
a boundary condition, our analysis establishes that the only way one can calculate the free-space 
radiation field of such sources is via the retarded solution for the potential. Finally, we discuss the 
applicability of these findings to pulsar observational data: the more distant a pulsar, the narrower 
and brighter its giant pulses should be. 
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I. INTRODUCTION 



Moving sources of electromagnetic radiation whose speeds exceed the speed of light in 
vacuo have already been generated in the laboratory UQ. These sources arise from sep- 
aration of charges: their superluminally moving distribution patterns are created by the 
coordinated motion of aggregates of subluminally moving particles. A polarization current 
density is, however, on the same footing as the current density of free charges in the Ampere- 
Maxwell equation, so that the propagating distribution patterns of such polarization currents 
radiate, as would any other moving sources of the electromagnetic field (5-^]. 

We have already shown, by means of an analysis based on the classical expression for 
the retarded potential [Eq. below], that the radiation field of a superluminally rotating 
extended source at a given observation point P arises almost exclusively from those of its 
volume elements that approach P, along the radiation direction, with the speed of light 
and zero acceleration at the retarded time js-10]. These elements comprise a filamentary 
part of the source whose radial and azimuthal widths become narrower (as 5r ~ Rp~ 2 and 
Sep ~ Rp~ 3 , respectively), the larger the distance Rp of the observer from the source, and 
whose length is of the order of the length scale l z of the source parallel to the axis of rotation 
[lo| . (r, ip, and z are the cylindrical polar coordinates of source points.) 

Once a source travels faster than its emitted waves, it can make more than one retarded 
contribution to the field that is observed at any given instant. This multivaluedness of the 
retarded time means that the wave fronts emitted by each of the contributing volume 

elements of the source possess an envelope, which in this case consists of a two-sheeted, 
tubelike surface whose sheets meet tangentially along a spiraling cusp curve (see Figs. 1 



and 4 of Ref . 



10|). For moderate superluminal speeds, the field inside the envelope receives 



contributions from three distinct values of the retarded time, while the field outside the 
envelope is influenced by only a single instant of emission time. Coherent superposition of 
the emitted waves on the envelope (where two of the contributing retarded times coalesce) 
and on its cusp (where all three of the contributing retarded times coalesce) results in not 
only a spatial, but also a temporal focusing of the waves: the contributions from emission 
over an extended period of retarded time reach an observer who is located on the cusp during 
a significantly shorter period of observation time. 

The field of each contributing volume element of the source is strongest, therefore, on 
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the cusp of the envelope of wave fronts that it emits. The bundle of cusps generated by 
the collection of contributing source elements (i.e., by the filamentary part of the source 
that approaches the observer with the speed of light and zero acceleration) constitute a 
radiation subbeam whose widths in the polar and azimuthal directions are of the order 
of S9p ~ -Rp" 1 and 5(pp ~ -Rp" 3 , respectively [l0(|. (_Rp, tpp and Op are the spherical 
polar coordinates of the observation point P.) The overall radiation beam generated by the 
source consists of a (necessarily incoherent ll|) superposition of such subbeams, a beam 
whose azimuthal width is the same as the azimuthal extent of the source and whose polar 
width arccos[c/(r<u;)] < \6p — ir/2\ < arccos[c/(r>u;)] is determined by the radial extent 
1 < r< < r < r> of the superluminal part of the source jsl, Q . (c is the speed of light in 
vacuo, uj is the angular frequency of rotation of the source, and f = ruj/c.) 

Since the cusps only represent the loci of points at which the emitted spherical waves 
interfere constructively (i.e., represent wave packets that are constantly dispersed and re- 
constructed out of other waves), the subbeams generated by a superluminal source need not 
be subject to diffraction as are conventional radiation beams. Nevertheless, they have a 
decreasing angular width only in the polar direction. Their azimuthal width b~(pp decreases 
as -Rp" 3 with distance because they receive contributions from an azimuthal extent 5tp of 
the source that likewise shrinks as Rp~ 3 . They would have had a constant azimuthal width 
had the azimuthal extent of the contributing part of the source been independent of Rp. On 
the other hand, the solid angle occupied by the cusps has a thickness 5zp in the direction 
parallel to the rotation axis that remains of the order of the height l z of the source distri- 
bution at all distances (see Fig. 2 of Ref. Q). Consequently, the polar width 56p of the 
particular subbeam that goes through the observation point decreases as -Rp -1 , instead of 
being independent of Rp [l Of ] . 

Because it is of a constant linear width, parallel to the rotation axis, an individual sub- 
beam subtends an area of the order of Rp, rather than Rp 2 . In order that the flux of 
energy remain the same across all cross sections of the subbeam, therefore, it is essential 
that the Poynting vector associated with this radiation correspondingly decay more slowly 
than that of a conventional, spherically decaying radiation: as -Rp -1 , rather than -Rp" 2 , 
within the bundle of cusps that emanate from the constituent volume elements of the source 
and extend into the far zone. This result, which also follows from the superposition of the 
Lienard-Wiechert fields of the constituent volume elements of a rotating superluminal source 
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@, Q] , has now been demonstrated experimentally . 

The narrowing of the individual subbeams with distance suggests that the absolute value 
of the gradient of the radiation field described here should increase with distance, in contrast 
to that of a conventional, diffracting radiation beam that decreases with distance. This is 
illustrated by a simple example. Imagine a rotating radiation beam with the amplitude 

A{R P ,cp P ,t P ) = A R P /2 exp[-( J R P ^ P ) 2 ], 

where -Rp stands for the scaled distance Rpu /c, </?p = ipp — utp is the azimuthal angle in 
the rotating frame, tp is the observation time, and Aq is a constant. This beam would be 
observed as a Gaussian pulse that has an azimuthal width of the order of Rp 3 and carries a 
constant flux of energy, 
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J A 2 # P 2 sin0pd0 P d^p = (2tv) 1/2 (c/uj) 2 A 2 

across any large sphere of radius Rp. The gradient of the amplitude of this pulse, 

dA/dipp = -2A R 7 p /2 (R P <pp) exp[-( J Rp^ P ) 2 ], 

increases in magnitude with distance as Rp 7 ^ 2 at the edges of the pulse. 

In this paper, we derive the azimuthal (equivalently, temporal) gradient (d/dtpp) of the 
radiation field that is generated by a physically viable, rotating superluminal source directly 
from the retarded potential, and show that the absolute value of this gradient does increase as 
R^ 2 within each subbeam. The spiky structure of the angular distribution of the emission 
from an accelerated superluminal source therefore follows not only from the geometry of 
the emitted cusps (geometrical optics) that was considered in Ref. [10|], but also from the 
calculation of the field distribution (physical optics) that is presented here. This result 
corroborates the earlier finding that the overall radiation beam consists of an incoherent 
superposition of sharply peaked subbeams that become narrower with distance from the 



source 
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There is, however, another, more significant implication. The retarded solution to the 
wave equation that governs the electromagnetic potential in the Lorenz gauge [Eq. (T5]) below] 
generally entails three terms: an integral over the retarded value of the electric current 
density, an integral over the boundary values of the potential and its gradient, and an 
integral over the initial values of the potential and its time derivative [see Eq. ([3]) below]. 



For a localized distribution of electric current, the integral over the retarded value of the 
source density is of the order of -Rp 1 in the far zone. If evaluated for a potential that is of 
this order of magnitude in the far zone (i.e., decays as -Rp 1 ), the integral over the boundary 
in this solution would also be of the order of -Rp 1 in the limit where the boundary tends to 
infinity. However, even potentials that satisfy the Lorenz condition are arbitrary to within 
a solution of the homogeneous wave equation, so that one can always use the gauge freedom 
in the choice of potential to set this boundary term identically equal to zero. 

In the case of the corresponding retarded solution of the wave equation for the electro- 
magnetic field [Eq. (JTj) below] , on the other hand, one no longer has the freedom offered by a 
gauge transformation to render the boundary term equal to zero. Nor does this term always 
decay faster than the source term, so that it could be neglected for a boundary that tends 
to infinity, as is commonly assumed in textbooks (e.g., page 246 of [ljj]) and the published 



literature 
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16(| . The boundary contribution to the retarded solution of the wave equation 



governing the field entails a surface integral over the boundary values of both the field and 
its gradient [see Eq. ([8]) below]. In the superluminal regime, where the gradient of the field 
increases as Rp 2 over a solid angle that decreases as -Rp -4 , this boundary contribution turns 
out to be of the order of _R p 1//2 (see Section HV|) . Not only is this not negligible relative to 
the contribution from the source term, which decays as .Rp 1 13Ml6l| . but the boundary term 
constitutes the dominant contribution toward the value of the radiation field in this case. 

Thus, if one ignores the boundary term in the retarded solution of the wave equation 
governing the field (as is done by Hannay 13l4l6l|). one would obtain a result, in the super- 
luminal regime, that contradicts what is obtained by calculating the field via the retarded 
potential 8KH|. However, the contradiction stems solely from having ignored a term in the 

A 1/2 

solution to the wave equation that is by a factor of the order of _R P greater than the term 
that is normally kept in this solution. The contradiction disappears once the neglected term 
is taken into account: the solutions to both the wave equation that governs the potential 
and the wave equation that governs the field predict that the field of a rotating superluminal 
source decays as _R p 1//2 as -Rp tends to infinity. 

From a physical point of view, however, what one obtains by including the boundary term 
in the retarded solution to the wave equation that governs the field is merely a mathematical 
identity; it is not a solution that could be used to calculate the field arising from a given 
source distribution in free space. Unless its boundary term happens to be negligibly smaller 
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than its source term, a condition that cannot be known a priori, the solution in question 
would require that one prescribe the field in the radiation zone (i.e., what one is seeking) as 
a boundary condition. The role played by the classical expression for the retarded potential 
in radiation theory is clearly much more fundamental than that played by the corresponding 
retarded solution of the wave equation governing the field. The only way to calculate the 
free-space radiation field of an accelerated superluminal source is to calculate the retarded 
potential and differentiate the resulting expression to find the field (see also We must 

again emphasize that this is an important contrast with subluminal sources. 

This paper is organized as follows: Section [III presents the retarded solutions to the 
initial-boundary value problems for the wave equations that govern the potential and the 
field. We provide a detailed mathematical derivation of the gradient of the radiation field 
that is generated by a rotating superluminal source in Section IHIt with a brief account of 
the required background material in Subsection IHI Al the formulation of the problem in 
Subsection IIIIBl the derivation of an integral representation of the gradient of the Green's 



function in Subsection IHI CI the regularization of the integral over the radial extent of the 
source in Subsection IHI Dl (and Appendix |A]), a description of contours of steepest descent 
in Subsection IIIIE| and the asymptotic evaluation of the gradient of the radiation field in 
Subsection IIIIFI Section [TV] evaluates the boundary term in the retarded solution to the 
wave equation governing the field, and we conclude in Section [V] 

II. BOUNDARY TERM IN THE SOLUTION TO THE WAVE EQUATION 

In the Lorenz gauge, the electromagnetic fields 

E = -V P A° - dA/d(ct P ), B = V P x A, (1) 

are given by a four-potential that satisfies the wave equation 

1 d 2 A> 1 4tt 

^~- 2 ^ = --f, A* = 0,-,3, (2) 
cr at 1 c 

where A°/c and j°/c are the electric potential and the charge density, and A 11 and for 
H = 1,2,3 are the components of the magnetic potential A and the current density j in a 



Cartesian coordinate system 12]. The solution to the initial-boundary value problem for 
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Eq. (|2J) is given by 



1 /' /: ' . f 1 



^(xp,tp)=-/ dt d i xfG + — dt dS- (GV A"- A"VG) 
cJo Jv ' 47r Jo is 



p 



d^a; #- G 



(3) 



Aire 2 Jv \ dt dt 

in which G is the Green's function and £ is the surface enclosing the volume V (see, e.g., 
page 893 of |lfi|). 

The potential that arises from a time-dependent localized source in unbounded space 
decays as i?p _1 when Rp ^> 1, so that for an arbitrary free-space potential the second term 
in Eq. ([3]) would be of the same order of magnitude (~ -Rp -1 ) as the first term in the limit 
that the boundary £ tends to infinity. However, even potentials that satisfy the Lorenz 
condition V • A + c~ 2 dA° /dt = are arbitrary to within a solution of the homogeneous 
wave equation: the gauge transformation 

A -> A + VA, A -> A - dA/dt (4) 



preserves the Lorenz condition if V 2 A — c~ 2 d 2 A/dt 2 = (see [l2j). One can always use 
this gauge freedom in the choice of the potential to render the boundary contribution (the 
second term) in Eq. ([3]) equal to zero, since this term, too, satisfies the homogenous wave 
equation. Under the null initial conditions A^\ t =o = {dA^ /dt) t =o = 0, assumed in this paper, 
the contribution from the third term in Eq. ([3]) is identically zero. 

In the absence of boundaries, the retarded Green's function has the form 

n 5(t P -t - R/c) 
G(x, t; x P , t P ) = lP (5) 

where 5 is the Dirac delta function and R is the magnitude of the separation R = xp — x 
between the observation point xp and the source point x. Irrespective of whether the 
radiation decays spherically or nonspherically, therefore, the potential A^ due to a localized 
source distribution, which is switched on at t = in an unbounded space, can be calculated 
from the first term in Eq. (J3j): 

A"(x P , t P ) = c' 1 J d 3 xdt j"(x, t)6(t P - t - R/c)/R, (6) 

i.e., from the classical expression for the retarded potential. Whatever the Green's function 
for the problem may be in the presence of boundaries, it would approach that in Eq. (jSJ) in 
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the limit where the boundaries tend to infinity, so that one can also use this potential to 
calculate the field on a boundary that lies at large distances from the source. 
Next, let us consider the wave equation that governs the magnetic field 

1 d 2 B 4vr . 

This may be obtained by simply taking the curl of the wave equation for the vector potenial 
[Eq. (j2J) for n = 1,2,3]. We write the solution to the initial-boundary value problem for 
Eq. (J7J), in analogy with Eq. (J3J), as 

£? fc (x P , t P ) =- [ P dt [ d 3 x (V x j) fc G + -j- [ P dt [ dS • (GVB k - B k VG) 

C Jo Jv 47r JO Jt. (o\ 



Anc 2 Jy V dt dt J t=o 

where k = 1,2,3 designate the components of B and V x j in a Cartesian coordinate system. 
Here, we no longer have the freedom, offered in the case of Eq. ([2D by a gauge transformation, 
to make the boundary term zero. 

Our task in this paper is to demonstrate that the boundary contribution in Eq. (jSJ) is, in 

" 1/2 

fact, by a factor of the order of R p larger than the source term of this equation in the far 
zone when the source is superluminal and accelerated. For this purpose, we need to know 
how the the gradient V-Bfc in the second term in Eq. (jBJ) decays in the far zone. We shall 
calculate, in the following section, the field B and its gradient directly from the classical 
expression for the retarded potential [Eq. (jSJ)], and use the resulting expressions to evaluate 
the second term in Eq. (jSj) for a boundary that lies in the far zone. 

III. GRADIENT OF THE RADIATION FIELD GENERATED BY A ROTATING 
SUPERLUMINAL SOURCE 

A. Background: The exact expression for the radiation field 

We base our analysis on the generic superluminal source distribution considered in 

PI r i \1 

Refs. [9] and [10], which has already been created in the laboratory [2j|. This source comprises 
a polarization current density j = dP / dt for which 

Pr,<p,z( r , f, z, t) = s r>tPtZ (r, z) cos(m<p) cos(fit), — n < (p < ir, (9) 
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with 



tp = cp — ut, 



(10) 



where P r ,ip,z are the components of the polarization P in a cylindrical coordinate system 
based on the axis of rotation, s(r, z) is an arbitrary vector that vanishes outside a finite 
region of the (r, z) space, and m is a positive integer. For a fixed value of t, the azimuthal 
dependence of the density ([8]) along each circle of radius r within the source is the same 
as that of a sinusoidal wave train, of wavelength l^rjm whose m cycles fit around the 
circumference of the circle smoothly As time elapses, this wave train both propagates around 
each circle with the velocity ru and oscillates in its amplitude with the frequency Q. This 
is a generic source: one can construct any distribution with a uniformly rotating pattern, 
P r ,<p,z( r , <P, z), by the superposition over m of terms of the form s rtip>z (r, z, m) cos{m(p). 

To find the retarded field that follows from Eq. (jSJ) for the source described in Eq. fl9]), we 
first calculated in Ref. the Lienard-Wiechert field of a circularly moving point source with 
a speed rut > c, i. e., a generalization of the synchrotron radiation to the superluminal regime. 
We then evaluated the integral representing the retarded field (rather than the retarded 
potential) of the extended source (jBJ) by superposing the fields generated by the constituent 
volume elements of this source, i.e., by using the generalization of the synchrotron field as 
the Green's function for the problem (see also 17]). In the superluminal regime, this Green's 
function has extended singularities that arise from the coherent superposition of the emitted 
waves on the envelope of wave fronts and its cusp. 

Inserting the expression for j = dP/dt from Eq. into Eq. (jHJ), and changing the 
variables of integration from (x, t) = (r, p, z, t) to (r, ip, z, 0), we found in Eq. (20) of Ref. ^ 
that the magnetic field B of the generated radiation is given by 

B = -|i(uVc) 2 £ M=(U± f v rdrd<pdzfiexp(-itpp) J2 3 j= i n j dG j/ d( P, (H) 

where fj,± = {0,/uj) ± m, 

m = s r cos#pe|| + Stp&x, U2 = — s^cos^pen + s r e±, U3 = — s z sin 9pe\\ , (12) 

and Gj (j = 1, 2, 3) are the functions resulting from the remaining integration with respect 
to <p: 



G 2 
G 3 



dip 



R 



■ exp(iQ<p / uj) 



cos(p — pp) 
sin(<£> — (p-p) 
1 



(13) 
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Here stands for ip — ipp with ipp = <^p — cutp, i? is 

.R = [(^p - z) 2 + r P 2 + r 2 - 2r P r cos(y? P -(p)]%, (14) 
the function g is defined by 

p = <p - ip P + R, (15) 

with R = Ru/c, Aip is the interval of azimuthal angle traversed by the source, and V is the 
volume occupied by the source in the (r, <p, z) space. The unit vector ey = (e z x h)/\e z x n| 
(which is parallel to the plane of rotation), e_i_ = nxeii, and the radiation direction n = R/i? 
together form an orthonormal triad (e 2 is the base vector associated with the coordinate 
z). The corresponding expression for the electric field in the limit Rp = |xp| — > oo, where 
fi ~ xp/|xp|, is given by E = fix B, as in any other radiation. 

A distinctive feature of the emission from a superluminal source is the multivaludeness of 



the retarded time 



tUio 



101 ]. At any given observation time, at least three distinct contributions, 
arising from three differing retarded times, are made toward the value of the radiation field 
by the part of the source that lies within the following volume of the (r, ip, z) space: 

A > 0, 0_ < <p < <f> + , (16) 

where 

A = (f|-l)(f 2 -l)-(i-z P ) 2 , (17) 



<p± = 2tt - arccos[(l =f Aa)/(ff P )] + R±, (18) 

and 

R± = [(z - zpf + f 2 + f 2 - 2(1 t A3)]3. (19) 

This volume is bounded by a two-sheeted surface, the so-called bifurcation surface, whose 
two sheets = 0±(r, z) meet tangentially along a cusp (see Figs. 3 and 4 of Ref. [9J). The 
strongest contributions are made by the source elements that lie close to the cusp curve 
A = 0, cf) = 4>±\a=o, where the two sheets of the bifurcation surface meet tangentially. For 
Rp 3> 1, the filamentary locus of these contributing source elements is essentially parallel 
to the rotation axis and has exceedingly narrow radial and azimuthal widths, of the orders 
of -Rp 2 and -Rp 3 , respectively (see Fig. 2 of Ref. H) 
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The asymptotic values of the Green's functions Gj close to the cusp curve of the bifur- 
cation surface (where A<C 1) are given by 

[ G/ n \x\ < 1 
Gj = I (20a) 

[or' \x\>h 

with 

Gj 111 ~ 2cx _2 (l — % 2 )~ 5 [pj cos(| arcsinx) — c\qj sin(| arcsinx)], (20b) 

and 

Gj OUt ~ Cx _2 (x 2 — l)~a [pj sinh(|arccosh|x|) + c\qjSgn{x) sinh(|arccosh|x|)], (20c) 
where 

X = 3(0-c 2 )/(2 Cl 3 ), (21) 

with 

d = - c 2 = §(</>+ + (22) 

and the symbol ~ denotes asymptotic approximation. The derivation of these asymptotic 



values, together with the exact expressions for the coefficients Pj(r, z) and qj(r,z) may be 
found in the Appendices of Refs. [8] and [9]. Here, we only need the following limiting values 
of these coefficients for Rp 1: 

Pi ~ 2^(w/c) J Rp 2 exp(i^ c /cj), (23) 

p 2 ~ -Rppi, p 3 ^ -P2, (24) 

and 

gi ~ 2^(cu/c).Rp 1 exp(ifi^ c /co>), (25) 



g 2 ~ -g 3 ~ -i(VL/uj)q u (26) 

where y? c ~ (pp + 3n/2 in this limit. Note that, in these expressions, G* n,out represent 
the different forms assumed by the Green's functions Gj inside and outside the bifurcation 
surface, i.e., for <p inside and outside the interval (</>_, </>+) respectively (see Fig. 6 of Ref. 9]]). 
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The above results show that as a source point (r, ip, z) in the vicinity of the cusp curve 
A = 0, (f> = 0±|a=o> approaches the bifurcation surface from inside, i.e., as x ^ 1- or 
X — > — 1+, Gj m and hence Gj diverge. However, as a source point approaches one of the 
sheets of the bifurcation surface from outside, Gj tends to a finite limit: 

G7 ut U ± = Cr\ =±1 * iPi ± 2c^)/(3 Cl 2 ); (27) 

for, the numerator of Gj OUt is also zero when \x\ = 1. The Green's function Gj is singular, 
in other words, only on the inner side of the bifurcation surface (see Fig. 6 of Ref. {9J). 



B. Formulation of the problem 

It turns out that none of the componenets of the gradient of B can be evaluated for the 
source distribution (Q without a lengthy calculation. However, we shall see in Section IIVI 
that the radial component of is of the same order of magnitude in the far zone as the 
azimuthal (or equivalently, temporal) component dB k /dipp of the gradient of B k (r P , <p P , z P ). 
Since this component of the field gradient is both algebraically simpler to calculate and more 
directly related to the observeable characteristics of the generated subbeams (Section HJ), it 
will be the only component that we shall here evaluate in detail. The relationship between 
the far-field values of this and the other components of the field gradient is not difficult to 
establish (Section HVl) . 

The component dH/d0p of the gradient of B may be calculated by differentiating the 
right-hand side of Eq. ( TTTT) under the integral sign and using the fact that dG/dip-p = 
—dG/dp. It follows from an argument identical to that given in Ref. 9J (in connection with 
calculating B itself) that the contribution (<9B/<9<£>p)a>o arising from the source elements in 
A > toward the value of dH/dtpp can be written as 

(0B/0£p) A > o = (dB/dif P ) m + (dB/d0 P ) out (28) 



with 



where 



(«9B/^ P r out = |i(^/c) 2 £ J =1 f A > r dr dz u, L>> ou \ (29a) 



(29b) 
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and 



out 



E 



+ 

7T — 0p J (f>+ 



TT-ifip 



d(p n exp(— ifup) (d 2 Gj/dip 



2\out 



(29c) 



Once it is integrated by parts, the integral in Eq. (129bj) in turn splits into three terms: 



(30) 



d0 exp(— ifiip)Gj 



of which the first two (integrated) terms are divergent [see Eq. (120bj) ]. Hadamard's finite part 
of Lj m and hence of (dB / d0p) m , here designated by the prefix J 7 , is obtained by discarding 
this divergent contribution toward the value of Lj m (see Refs. {9] and 



E 



<f>+ 



exp(—ifj,if)Gj m . 



(31) 



Note that the singularity of the kernel of this integral, i.e., the singularity of G/ n , is like 
that of \<p± — y3| - 5 and so is integrable. 

The boundary contributions from = <f>± that result from the integration of the right- 
hand side of Eq. (129c[) by parts are well-defined automatically: 



U 



{/iexp(-i^) UdGj/dfi)™* + ifiGj 



4>+ 



|U=£t± 



+ 



<t>- rir-tpp . 

+ / )d0/i 3 exp(— ifi0)Gj 

■TV — ipp J (f>+ 



(32) 



since (dGj/d0) ont (like Gj OUt ) tends to a finite limit as the bifurcation surface is approached 
from outside (see Subsection IIII Cj) . In deriving Eq. ( 132|) . we have made use of the fact 
that (dGj /dipY^l^-^p equals (dGj/dcp) 0111 ]^^-^ when <p± ^ ±n — (pp. The integral 
representing Lj OU , in other words, is finite by itself and needs no regularization. 

If we now insert F{L/ n } and L/ Ut from Eqs. §H) and ([32]) in Eq. fl29a| and combine 
(dB/d(pp) m and (dB/d0p) out , we arrive at an expression for the Hadamard finite part of 
a>o which entails both a volume and a surface integral: 

F{(dB/d0 P ) A > o } = (dB/dcppf + («9B/<^ P ) ns . (33) 

The volume integral 

{dB/d0 P ) s = - \i(u:/c) 2 En=n ± ^ Ia>o r dr dz /- w d< £ exp(-w) 

(34) 



x 
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las the same form as the familiar integral representation of the field of a subluminal source 
I2I and decays spherically (as i?p _1 for Rp ^> 1). 
The surface integral 

(dB/d£rY s = -\^/cf Y.U L>o ' dr dz ( 35 ) 
stems from the boundary contribution 

V dgC = Yl /" ex P(-W) 

in Eq. fl32|) . It is this contribution that turns out to increase, rather than decay, in the limit 
i?P — ► 00. To see this, we need to know the values of (dGj/dip) out at <p = <fr±, in addition to 
those of Gj° nt \<j >=l f )± which are given in Eq. ( 1271) . 



(dGj/d<p) out + ifuGj 



(36) 



C. Azimuthal (or temporal) gradient of the Green's function 



The Green's function (1131) depends on (p only through the variable <fi which appears in 
the argument of the Dirac delta function, so that the differentiation of Eq. ( jTBl) with respect 
to <p simply yields 



dGj/d<p 



diphj(if)8'(g - 0), 



(37) 



where 5' stands for the derivative of the delta function with respect to its argument, and 



hi 
h 2 
h 3 



exp(iQ{p/ui) 
R 



cos(<p — y?p) 
sin(y9 — ip P ) 
1 



(38) 



Integrating the right-hand side of Eq. (!37j) by parts, we obtain 



d_G± 

dip 



A(p og/ dip dip 



dip — 



hj(<p) 



(39) 



-dg/dip. 

when the source trajectory intersects the bifurcation surface of the observation point (i.e., the 
argument of the delta function vanishes within Aip). A uniform asymptotic approximation 
to this inte gral , for small A, can be found by the method of Chester et al. in the time 



domain 



20 
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Where it is analytic (i.e., for all x 7^ xp), the function g(p>) transforms to the cubic 
function 

g{p>) = |z/ 3 - Cl 2 z/ + c 2 , (40) 

where v is a new variable of integration replacing p and the coefficients c\ and c 2 [defined in 
Eq. (1221) ] are such that the values of the two functions on opposite sides of Eq. fj4T)]) coincide 
at their extrema. Insertion of Eq. (140]) and its derivative, 

dg v 2 — c 2 



in Eq. (J39J) results in 



dG 

dp> 



' - 1 dv 

Au 



dip dip jdv 



F'a 



[y 2 — Ci 2 ) 2 V 2 — Ci 2 



5(|z/ 3 - c 2 v + c 2 - 



where 



] -\dv) dp 2 i! 



(41) 



(42a) 



(42b) 



and Az/ is the image of Ap under transformation (j40]). 

As in the evaluation of Gj in Refs. |8[ and [9] , the leading term in the asymptotic expansion 
of the integral (I42ap for small ci, which corresponds to small A [see Eq. fl53l) below], can 
now be obtained by replacing the functions F and F 1 in its integrand with Pj + Qji/ and 
P'j + Q'jV, respectively, and extending its range Az/ to (— 00, 00): 



dp 



00 

~ ' dz/ 

00 



P. + Q,v P'j + Q'jU 



(v 2 - Cl 2 ) 2 



+ 



V 2 — Ci 2 



#Qz/ 3 - c 2 v + c 2 



where 



-c\/i 



(43a) 



(43b) 



Qj 2*^1 (^j|i / =ci Fj\u=— ci); 



(43c) 



and 



^ j ~~ 2 (-^ i I ^= c i "I" F j I i/=-ci ) > 



2^1 (^il^=ci F j\v=— ci)- 



(43d) 



(43e) 
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Note that the extrema 

ip± = 2-k - arccos[(l =F A^)/(ff P )] (44) 

of the function g(ip) transform into v = =)=ci, respectively. 

The derivatives dip / dv\ v= ± Cl that appear in the definitions of the coefficients 
(Pj,Qj,P'j,Q'j) are indeterminate. Their values must be found by repeated differentia- 
tion of Eqs. ffl~5]) and (40) with respect to u: 

(dg/d(p)(d(p/dv) = v 2 - ci 2 , (45a) 



(d 2 g/d^ 2 )(d V /du) 2 + (dg/d<p)(d 2 <p/dv 2 



2za 



(45b) 



etc., and the evaluation of the resulting relations at v = ±c\. This procedure, which amounts 
to applying l'Hopital's rule, yields 



d(p/dv\ u= ± ci = (2c iJ R T )2/A 4 . 



Using d 2 g/d(p 2 



T& 1/2 /R± and Eq. (g6]), we find from Eq. fl42b|) that 



I v=±c\ 



±2c x /, 



li/=±ci' 



(46) 



(47) 



in which fj = (d(p/di>)hj are the functions earlier encountered in the evaluation of Gj in 
Refs. |8| and [9]. Hence, Pj = 1c 2 and Qj = 2pj, where pj and qj are precisely the same 
as the coefficients in Eqs. (T2"U|) that are approximated in Eqs. (T2"51) - (|2"6"I) (see Ref. [9]). 

We now need to evaluate dGj/dip only outside the bifurcation surface, i.e., for \x\ > 1 
[see Eqs. f l2Tj) and fl36l) ]. In this region, the argument of the S function in Eq. f)43ap has a 
single zero at 

v = v* = 2cisgn(x) cosh(|arccosh|x|), |x| > 1, (48) 
(see Appendix A of Ref. [§]). Outside the bifurcation surface, therefore, Eq. (143 al) yields 



dGj\^t 



80 



\h> 2 — Ci 2 | 



Pj + QjU P'j + Q'ju 



y - Cl 2 ) 2 



V 2 — C\ 2 



(49) 



Keeping only the first term in this expresssion, which is dominant when c\ <C 1, we obtain 
/dGj\ out 2sinh 3 (iarccosh|x|) r ,-, -, 

\~W ) ~ Ci 5 (x 2 ~ 1) 3/2 ^ + 2PjSg ^ COsh (3 arCCOSn ^l)J ' ( 50 ) 

in which Pj(r, z) and qj(r, z) assume the values given in Eqs. (T23T) f[26|) when R P 1. 
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Evaluation of the right-hand side of Eq. (|50|) at (j) = <f>± now yields the following term 
that appears in the expression for Lj edge : 



expi-iwKdGj/dft^r ~ - (|) 3 c 1 - 5 exp[-i / i(^ P + c 2 )] 



2 \ 3 _5 

■,-A.^| — j/M ^ | > -r ( o ; | 

(51) 



x [pj cos(|/iCi 3 ) — |ici9j sin(|/xci 3 

The asymptotic expansions of G° nt and (dGj/dip) ont given in Eqs. f!20cj) and (!50|) are for 
small C\. To be consistent, we must likewise replace the expression that is found by inserting 
Gj OUt and (dGj/d0) out in Eq. ( l36l) with the leading term in its expansion in powers of C\. 
The result is 

L .edge „ _2 1 /3 ( | ) 3 jRp -l^. Ci -5 exp(i ^ c/L , ) ^ eX p [-i/i^P + </>_)], (52) 

in which pj = (Rp 1 — 1 +1) and y? c ~ y?p + 37r/2. 

The far-field value of c\ close to the cusp curve of the bifurcation surface (where A = 0) 
is given by 

ci ~ 2 _ ^p X A3 (53) 
[see Eq. fl22l)]. Inserting Eq. (J52"|) in Eq. fl35|) and using Eq. fl53|) . we finally arrive at 

(9B/^ P ) ns ^2(|) 3 J R P exp[i(^ c /^ - tt/2)] ^ =w/i exp(-# P ) 

3 " _5 (54) 

fdfd^A 2Uj exp(— i/i0_). 

>o 



3 

Eft / 

,=i 



X 

r- 

As in Ref. [10|, the integral over f in this expression may be evaluated by contour integration. 
Since the singularity of its integrand at A = is not integrable, however, the contour integral 
that passes through this singularity needs to be in addition regularized. 



D. Regularization of the integral over the radial extent of the source 

The kernel of the integral in Eq. (|54[) has the same phase but a different amplitude 
compared to that of the integral encountered in Eq. (19) of Ref. [l(3|. Hence, the asymp- 
totic evaluation of integral (1341 entails the same techniques as those used before, but the 

n 

regularization of this integral requires an extension of the procedure followed in Ref. [101 ] . 

The function 0_(f, z) in the phase of the integrand in Eq. (|54|) is stationary as a function 
of f at 

f = f c (z) = {|(f| + 1) - [\{fl - l) 2 - (z - ip) 2 ] 1 /*}!/ 2 . (55) 
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When the observer is located in the far zone, this one-dimensional locus of stationary points 
coincides with the locus, 

r = rs=[l + {z-z?) 2 /{rl-l)] l l\ (56) 

of source points that approach the observer with the speed of light and zero acceleration 
at the retarded time, i.e., with the projection A = of the cusp curve of the bifurcation 
surface onto the (r, z) plane (see Fig. 4 of Ref. n]]). For Rp ^> 1, the separation fc — fs 
vanishes as Rp 2 [see Eq. (|69|) below] and both fc and fs assume the value csc#p. 
It follows from Eq. ffl8l) that at the stationary point f = fc, 

4>-\? =fc = 4> c = Rc + <Pc - <Pp, (57) 

d(p-/df\r=r c = 0, and 

d 2 <p./df% =fc = a = -R c l [{fl - l){fl - I)" 1 - 2], (58) 

where 

^Pc = V^p + 2vr — arccos(f c /f P ) and R c = f c (fp — f c )i . (59) 

For observation points of interest to us (the observation points located outside the plane of 
rotation, #p 7^ 7r/2, in the far zone, Rp ^> 1), the parameter a in Eq. (158]) has a value whose 
magnitude increases with increasing Rp: 

a~ -Rpsm 4 6psec 2 6p. (60) 

In other words, the phase function 0_ is more sharply peaked at its maximum, the farther the 
observation point is from the source. This property of the phase function <p_ distinguishes the 
asymptotic analysis that will be presented in this section from those commonly encountered 
in radiation theory. What turns out to play the role of a large parameter in this asymptotic 
expansion is distance (-Rp), not frequency (fi±u). 

The first step in the asymptotic analysis of the integral in Eq. (|5"1"1) is to introduce a 
change of variable £ = £(f , z) that replaces the original phase 0_ of the integrand by as 



simple a polynomial as possible j22|. This transformation should be one-to-one and should 
preserve the number and nature of the stationary points of the phase. Since 0_ has a single 
isolated stationary point at f = fc(z), it can be cast into a canonical form by means of the 
following transformation: 

<Mf, z) = <Mz) + ia(z)£ 2 , (61) 
19 



in which a is the coefficient given in Eqs. (1581) and (1601) . 
The integral in Eq. (l54"j) can thus be written as 

/ fdfd^A" 5/2 u i exp(-i / u0_) = / dM£ F(£, z) exp(ia£ 2 ), (62) 

in which 

5) = fA- 5 / 2 Uj (9f/^) exp(-i/i0 c ), (63) 

with 

df/di = <fiL(f 2 - 1 - A 1/2 )~\ (64) 

and a ee — fia/2. The stationary point f = fc and the boundary point f = f$ respectively 
map onto £ = and 

£ = & ee -[2^* - </> c )] 1/2 , (65) 

where 

<f) S ee 0_|, = , s = 2tt - arccos[l/(f 5 f P )] + (f|f 2 - I) 1 ' 2 . (66) 

The upper limit of integration in Eq. ( 1621) is determined by the image of the support of the 
source density (s in Uj) under the transformation (ISTj) . 

By substituting the value of r c from Eq. (|55|) in Eq. (flTI) . we find that A 1 / 2 = f 2 — 1 at C. 
Thus, the Jacobian df/dl; of the above transformation is indeterminate at £ = 0. Its value 
at this critical point must be found by repeated differentiation of Eq. (I6TT) with respect to £, 

(d(f)„/dr)(df/dO = a£ } (67) 



(<9 2 0_/<9f 2 )(<9f/<9O 2 + (<90_/<9f)(<9 2 f/<9£ 2 ) = a, (68) 

and the evaluation of the resulting relation ( |68l) at f = fc with the aid of Eq. ( 1581 . This 
procedure, which amounts to applying l'Hopital's rule, yields df/d£\^ = o = 1: a result we 
could have anticipated from the coincidence of transformation (l6Tj) with the Taylor expansion 
of 0_ about f = fc to within the leading order. Correspondingly, the amplitude F(^) in 
Eq. (J63l) has the value rc(r 2 — l)~ 5 \ij\r=r c exp(— ifi(pc) at the critical point C. 

To an observer in the far field (i?p 3> 1), the phase of the integrand on the right-hand 
side of Eq. (l62jl is rapidly oscillating, irrespective of how low the harmonic numbers /i± (i.e., 
the radiation frequencies /i±cu) may be. The leading contribution to the asymptotic value 
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of integral (162]) from the stationary point £ = can therefore be determined by the method 
of stationary phase. However, in the limit Rp — > oo, £$ reduces to 

£ s ~ -3~ 1/2 cos 4 6 P esc 5 # Pj R p 2 , (69) 

so that the stationary point £ = is separated from the boundary point £ = £s by an 
interval of the order of -Rp 2 only. We therefore need to employ a technique for the asymptotic 
analysis of integral ( 162]) that is capable of handling the contributions from both fc and f s 
simultaneously. 



E. Contours of steepest descent 



The technique we shall employ for this purpose is the method of steepest descents 22]. 
We regard the variable of integration in 

J(z) = r d£F(£,i)exp(k*£ 2 ) (70) 

Ha 

as complex, i.e., write £ = u + if , and invoke Cauchy's integral theorem to deform the 
original path of integration into the contours of steepest descent that pass through each of 
the critical points £ = £s, £ = and £ = £>. Here, we have introduced the real variable £>(i) 
to designate the image of f> under transformation (1ST]) , i.e., the boundary of the support 
of the source term Uj that appears in the amplitude F(£, z). We shall only treat the case in 
which /j, (and hence a) is positive; J(z) for negative fi can then be obtained by taking the 
complex conjugate of the derived expression and replacing 0c with —0c [see Eq. ( 163]) ]. 

The path of steepest descent through the stationary point C at which £ = is given, 
according to 

i£ 2 = -2uv + i(u 2 -v 2 ), (71) 
by u = v when a is positive. If we designate this path by C\ (see Fig. [Q, then 

f roc 

/ d£F(£, z) exp(ia£ 2 ) = (1 + i) / dvF\ (={1+i)v exp(-2av 2 ) 
JCi J-00 (72) 

~ (2tt//i) 1/2 exp ]-i(/i0 c - tt/4)] uj \ c sin 7 9 P \ sec 9 P \ 9 R~ 1/2 , 

for Rp ^> 1. Here, we have obtained the leading term in the asymptotic expansion of the 
above integral for large Rp by approximating F|^=(i+i) v by its value at C, where v — 0, and 
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using Eqs. (155]) and (!60|) to replace fc and a by their values in the far zone. Note that 

* —1/2 

the next term in this asymptotic expansion is smaller by a factor of order R p than this 
leading term. 



V 




FIG. 1: The integration contours in the complex plane £ = u + iv. The critical point C lies at the 
origin and us and u > are the images under transformation (|6ip of the radial boundaries f = rs(i:) 
and f = r>(z) of the part of the source that lies within A > 0. The contours C±, C*2, and C3 
are the paths of steepest descent through the stationary point C and through the lower and upper 
boundaries of the integration domain, respectively. 

The path of steepest descent through the boundary point S, at which u = us = £s and 
v = [see Eqs. fl65|) and fl69l) ]. is given by u = — (v 2 + U5 2 ) 1 ' 2 , i.e., by the contour designated 
as C2 in Fig. [TJ The real part of 

i£ 2 | C2 = 2^ 2 W) 1/2 W, (73) 

is a monotonic function of v and so can be used as a curve parameter for contour C2 in place 
of v . If we let 2v {y 2 + u^ 2 ) 1 / 2 = — r 2 , then it follows from 

£k = ~K 2 + ir 2 ) 1/2 (74) 

that 

/ d£F(£, 5) exp(i< 2 ) = exp[i(cm s 2 - tt/2)] 

JC2 ,00 (75) 

x J Q drr( M5 2 + ir 2 )- 1 / 2 F| 5= _ (us2+ . r2)1/2 exp(-ar 2 ). 
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The function F\c 2 in this expression has to be determined by inverting the following version 
of the original transformation (IBTj) : 

<Mf,z)-<Mz) = iiar 2 . (76) 

Here, we have used Eqs. (I65I) and (IT4I) to rewrite Eq. (161]) in terms of r. 

Since the dominant contribution toward the asymptotic value of the above integral for 
Rp ^> 1 comes from the vicinity of the boundary point S, where r = 0, the required 
inversion of transformation ( 1761) can be carried out by means of a Taylor expansion of the 
phase function 0_ in powers of r (see Appendix|A|). We find in Appendix |A]that the resulting 
expression for F(t)\c 2 diverges at r = as r~ 4 . Therefore, as in the case of the integral over 
if in Eq. ( JTTT) . we must regard the divergent integral in Eq. ( 1751) as a generalized function 



that equals its Hadamard's finite part (see, e.g., Ref. 13]). The procedure for finding the 
Hadamard finite part of this integral, though lengthy, is straightforward and results in 

f[J^ d£F(£,£)exp(ia£ 2 )} =(35/6 4 )(27r/ yU ) 1 / 2 exp[-i(^ s - 3vr/4)] 



c 2 > (77) 

.j | sin 7 6p | sec Op^Rp 1 ^ 2 



in the limit i?p>l (see Appendix lAl) . 

The path of steepest descent through the boundary point £ = £>, at which u = u > , v — 0, 
is given by k = (t> 2 + u>) 1//2 , i.e., by the contour designated as C 3 in Fig. [H The real part 
of the exponent 

i£ 2 |c 3 = -2v(v 2 + U> 2 ) 1/2 + i U> 2 (78) 

is again a monotonic function of v and so can be used to parametrize contour C3 in place of 
v. If we let 2v(v 2 + m> 2 ) 1 / 2 = %■> then it follows from 

e| C3 = ( M> 2 + i X ) 1/2 (79) 

that 

/ d£F(£, z) exp(i< 2 ) = \ exp[i(aM> 2 - vr/2)] 

dx (m> 2 + ix)" 1/2 ^| €=(u>2+ix) i/2 exp(-ax)- 



(80) 



The asymptotic value of this integral for Rp 3> 1 receives its dominant contribution from 
X = 0. Because the function F\ Cs is regular, on the other hand, its value at x = can be 
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found by simply evaluating the expression in Eq. (IBBl at r — f>. The result, for Rp — > oo, 
is 

F\c 3 , x =o — r 2 _Rp 4 sin 4 # P sec 2 #p(f 2 sin 2 # P — l) _3 u : ,| ?= ^ > exp(— i/i0 c )-u > (81) 
[see Eqs. (fTTl) and (I60l) ]. This in conjunction with Watson's lemma therefore implies that 

/ d£F(£, z) exp(i< 2 ) ~f 2 (f 2 sin 2 6 P - l)" 3 ^^ 

Jc 3 (82) 

x exp[-i(0_|r=r> + 7r/2)]/i~ 1 ^p 5 , 

to within the leading order in Rp . 

F. Asymptotic value of the gradient of the field for large distances 

The integral in Eq. (ITUI) equals the sum of the three contour integrals in Eqs. (!72|) . (!77j) . 
and (1821) ; the contributions of the contours that connect C\ and C*2, and Ci and C3, at 
infinity (see Fig. [1]) are exponentially small compared to those of Ci, C2, and C3 themselves. 
On the other hand, the leading term in the asymptotic value of the integral over C3 decreases 
(with increasing Rp) much faster than those in the asymptotic values of the integrals over 
C\ and Ci'- the integral over C3 decays as -Rp 5 , while the integrals over C\ and C2 decay 

^ 1/2 _____ 

as R p . According to Eqs. (154|) . (1621) . fl72l) and (1771) . the leading term in the asymptotic 
expansion of the contribution (<9B/<9<y3p) ns , for large Rp, is therefore given by 



(dB/d<p P ) ns -±(35 - 6 4 i)i? P /2 sin 7 flp| sec 6 P \ 9 exp{i[(n/u)<p c + tt/4]} 

(27r|/i|) 1/2 sgn(/i)exp(ifsgn/i) 



3 7 

x 



3 

X 

3 



(83) 



3 »oo 

/ diexp[-i^(^c + ^p)]uj|o, 
■ =1 7 -00 



in which /x± can also be negative (see the first paragraph of Section nil El) . 

The remaining z integration in the above expression for (<9B / '<9y3p) ns amounts to a Fourier 
decomposition of the source densities s r>tP)Z \c with respect to z. Using Eqs. (j~T) - (j~?l) to 



replace 4>c in Eq. ( 1831) by its far-field value 

(pc^Rp-z cos 6 P + 3vr/2, (84) 
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and using Eq. (fl2l) to write out uj in terms of s r>tPiZ , we find that 

{dB/d0 P ) ns ~^(35 - 6 4 i)i? P /2 sin 7 # P | sec# P | 9 

x exp{i[vr/4 + (Q/w)(<pp + 3vr/2)]} 

(85) 



J2 (27r| /U |) 1 / 2 sgn(/i) exp{i[f sgn(/i) - p(R P + <p P + 3vr/2)]} 



x 



x [(sp cos #p — s 2 sin ^p)ey — s r e±] , 
where s rj(/Pj2 stand for the following Fourier transforms of s riip>z \c with respect to z: 

/CO 
dz s r)V?)Z (f, z) | f=cscep exp(i^5 cos Op). (86) 
-oo 

Being the contribution from the source elements that approach the observer with the speed 
of light and zero acceleration at the retarded time, this expression is valid only at those 
polar angles 9 P within the interval arccos(l/f < ) < \9p — n/2\ < arccos(l/f > ) for which 
s r>iPjZ \f =CSC p is nonzero, i.e., at those observation points (outside the plane of rotation) the 
cusp curve of whose bifurcation surface intersects the source distribution. At these polar 
angles, the above expression for (<9B/<9y3p) ns constitutes the dominant contribution toward 
the gradient <9B/<9y3 P of the magnetic field of the radiation (see Subsection IIII Bp . 



IV. EVALUATION OF THE BOUNDARY TERM IN THE RETARDED SOLU- 
TION TO THE WAVE EQUATION GOVERNING THE FIELD 

Let the boundary E in the second term of Eq. (jSj) be a large sphere enclosing the source. 
The element dS of area for this boundary then has the form p 2 sin 6d6dipe p , where (p,<p,0) 
are the spherical polar coordinates in the space of source points, i.e., are related to the 
cylindrical polar coordinates (r, (p,z) we have been using by 

p = (r 2 + z 2 ) 1/2 , = arctan(r/^), (87) 

and e p is a unit vector in the direction of increasing p. Inserting this in the integrand of the 
boundary contribution in Eq. (jSj), we obtain 

B boundary = p 2 J dt J dip d6 sin (GdB/dRp \ Rp=p - BdG/dp) , (88) 

since (e p • V)B = d~B/dp = dTS/dRp\n p=p . We will be identifying the magnetic field B and 
its gradient on the boundary S with those of the radiation field that arises from source (Q. 
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These terms, which act as densities of two-dimensional sources in Eq. (|88|) . both have rigidly 
rotating distribution patterns, i.e., are functions of t in the combination (p = tp — ut only 
[see Eq. (183])]- 

Once the free-space Green's function fl5]) is inserted in Eq. ( I8"5j) . we can therefore cast it 
in the same form as Eq. (FTT]) by changing the integration variable t to ip, resulting in 

B boundary = p 2 J^dO sm6(G b dB/dR P \ Rp=p -BdG b /dp), (89) 

with 

G b = J dpR-^ig-cj)), (90) 

where p = pu/c, and g and are the same functions as those appearing in Eqs. ffT3l - ffT5l) . 
Equation (!90|) implies that 

~ R\dg, 

= \R + pRp sin 6* sin# P sm(ipj — ip P ) | _1 

where are the solutions of the transcendental equation g(cp) = [see Eq. (fl5!)]. For p ^> 1, 
the number of retarded positions <pj of the rapidly rotating distribution patterns of B|s and 
dH/dRp\£ that contribute toward the value of G b can be appreciably larger than three (see 
Ref. fl). 

The expression in Eq. fllip for the magnetic field B depends on Rp through dGjjd(p only, 
so that 

dB/dR P = - \\{yjjcf E M =^ ± fv rdrd<pdzfiexp(-i/j4>) 

3 (92) 
x 22 Uj&Gj/dRpdip. 

3=1 



(91) 



with 

d 2 G j /dR P dcp = -{u/c) [ dtphji^R^ld'^g - <p) - R-^'ig - <p)]dR/dR P , (93) 

where /ij are the functions defined in Eq. (I3"5j) . and a prime denotes differentiation of the 
delta function with respect to its argument [see Eq. (1371) ]. 

It follows from a comparison of the calculations described in Subsections IIII Al and IIII CI 
that the contribution of the second term on the right-hand side of Eq. (1931) toward the value 
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of dH/dRp is by a factor of the order of -Rp 4 smaller than that of the first term. Ignoring 
this small term, we obtain an expression for —d 2 Gj/dRpd<p that differs from the expression 
for d 2 Gj/d(p 2 only by the factor 

dRj dRp = FT 1 [Rp — z cos 6p — r sin 6p cos(<£> — <pp)\ (94) 

[see Eq. (1141) ]. which reduces to 1 in the limit Rp ^> 1. Correspondingly, the leading contri- 
bution, (d~B/dRp) ns , toward the value of dH/dRp is given by an expression identical to that 
in Eq. (IHo'j) for (<9B/cfy2 P ) ns , except that it is multiplied by —1 (see Section HTT|) . The absolute 
value of dH/dRp is therefore of the same order of magnitude as that of d 2 Gj/dRpd(p, and 



so increases as Rp 2 (see Subsection IIIIFI) . 



We already know that the radiation subbeams that are generated by the superluminal 



source ([9]) have the widths 50 ~ -Rp 1 and 5<p ~ -Rp 3 (see Ref. 10]). In the limit where the 
values of p and -Rp (i.e., the positions of the boundary and the observer) tend to infinity 
independently of each other, the Green's function G b reduces to 

Gb ~ [pRpsin8sin8p\sin(ifj — (pp)\j . (95) 

tp=tpj 

Hence, the absolute values of Gb and dG^/dp diminish with distance as (pRp)' 1 and p~ 2 Rp l , 
respectively. Since \dB/dRp\ increases as -Rp^ 2 while |B| decreases as Rp 1 ^ 2 , this means that, 
of the two terms inside the parantheses in Eq. (l89il , the second is negligibly smaller than the 
first. Inserting the orders of magnitude of the remaining factors in Eq. (I59"|) . in the order in 
which they appear in the first term of this equation, we obtain 

|B boundary | ~ p 2 x p~ 3 x p- 1 x (pRpy 1 x p 7/2 

( 96 ) 

~ P 1/2 Rp\ 

Thus, the absolute value of the boundary contribution toward the value of the field decays 

~ —1/2 

as Rp when the radius p of the spherical boundary S and the coordinate Rp of the 
observation point P are both large and of the same order of magnitude. 



V. CONCLUDING REMARKS 



The unaviodably lengthy calculation 
lends support to the conclusions of Ref. 



;hat we have presented in Sections [TTT] and [IV] both 
lo| , on the morphology of the radiation beam that 
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is generated by a polarization current with a superluminally rotating distribution pattern, 
and clarifies a fundamental issue concerning the method of calculating the radiation field 
of such a polarization current, which has been the source of a long-standing controversy 



13Nl71|. This calculation establishes 



(i) that the absolute values of both the radial component <9B / dRp and the azimuthal or 
temporal component <9B/cfy3p of the gradient of the radiation field that is generated 
by the superluminal source distribution (Q in the far zone are of the order of Rp 2 at 



any observation point within the overall radiation beam arising from this source 



23|, 



(ii) that the angular distribution of the emitted field contains sharply focused structures, 
i.e., that the overall radiation beam is composed of an incoherent superposition of 
rapidly narrowing subb earns jll]. 

(iii) that the boundary contribution toward the the solution of the wave equation governing 
the field decays as R p as the boundary tends to infinity, i.e., that the second term 
in Eq. (jSJ) is by a factor of the order of R p greater than the first term in this equation 
for a £ that lies in Rp ^> 1, 

(iv) that the discrepancy between the predictions of Eq. ([2]) and Eq. ([7]) disappears once 
one includes the boundary term that is normally neglected in solution ([8]), and 

(v) that Hannay's erroneous contention that the field of a rotating superluminal source 



should diminish as -Rp -1 , as does a conventional radiation field 13Ml6l|. stems from 
his having neglected the boundary term in the solution [Eq. ([8])] to the wave equation 
governing the field [Eq. (J7J)]. 

The sharply focused radiation pulses encountered in the present analysis are in fact 
observed in astronomical objects that are thought to contain superluminal sources. The 



radio emission received from pulsars is composed (often entirely 24|) of a collection of so- 



called giant pulses whose widths are as narrow as 1 ns |25] and whose brightness temperatures 



are as high as 10 39 K 2fjj]. Hankins et al 25[] note the puzzling brightness of these pulses: 



The plasma structures responsible for these emissions must be smaller than one 
meter in size, making them by far the smallest objects ever detected and resolved 
outside the Solar System, and the brightest transient radio sources in the sky. 
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The small size of the emitting structures reflects, in the present context, the narrowing (as 
Rp~ 2 and Rp~ 3 , respectively) of the radial and azimuthal dimensions of the filamentary part 
of the source that approaches the observer at P with the speed of light and zero acceleration 
at the retarded time n]]. This, together with the nonspherical decay of the individual 
subbeams generated by such filaments (as i?p~ 1//2 instead of i?p _1 ), easily accounts for the 
observationally inferred values of the brightness temperature of the giant pulses. 

The azimuthal (or temporal) gradient of the intensity of these pulses often appears in- 
finitely sharp at either their leading or trailing edges (see Fig. 1 of Ref. 3]). Correspond- 
ingly, the emission mechanism discussed in this paper sets no upper limit on the gradient 
d/d0p of the radiation field (i.e., on the sharpness of the leading or trailing edge of the 
pulse), if the length scale of spatial or temporal variations of its source are comparable with 
Rp 3 . According to the superluminal model of pulsars [^tJ (to which the present findings 
apply), the more distant a pulsar, the narrower and brighter its giant pulses should be. 
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APPENDIX A: HADAMARD'S FINITE PART OF THE DIVERGENT CONTRI- 
BUTION FROM THE INTEGRAL OVER THE CONTOUR C 2 

Our first task in this Appendix is to make the dependence of the integrand of integral 
f!75|) on the integration variable r explicit. This entails (i) inverting Eq. (1761) in the vicinity 
of the critical point r = to obtain f as a function of r for a fixed value of z, and (ii) 
expanding the function F\c 2 that appears in the integrand of integral (ITo'j) in powers of r. 
Next, we calculate the Hadamard finite part of the resulting integral (whose integrand turns 
out to diverge as r~ 4 at r = 0) by following the standard procedure used in the literature 
on generalized functions 
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Because it contains the factor 

A 1/2 = (f 2 - l) 1/2 (f 2 - r|) 1/2 , (Al) 

the function <p-(f, z) is not analytic at r — r$ [see Eqs. ( 1171) . (|18p and (|56j) ]. If, however, we 
eliminate f in 0_ in favour of 

n = {?-p s yi*, (A2) 

then the resulting function 

4>-{ V , z) ={r 2 P (r 2 s + r, 2 ) - [1 + (f 2 - l)Wri\ 2 }W 

+ 2tt - aiccos{fp\fl + r/ 2 )" 1/2 [l + (f 2 , - l) 1/2 r/]} 
can be expanded into a Taylor series about rj = to obtain 

0_ =0 5 + li^- 1 C os 2 6 P r] 2 - | sin 3 fl P r/ 3 + i^; 3 cos 2 ^ P (5 sin 2 6 P - l)r/ 4 
+ | sin 5 6> P ?7 5 H . 



(A3) 



(A4) 



Here, the coefficients in this series are approximated by their dominant values for Rp ^> 1, 
and the coordinate fg that appears in them is replaced by its value esc #p at the radius from 
which the main contributions toward the field in the far zone arise [see Eq. ( 156]) ]. Equation 
( 1A4[ ). in conjunction with Eq. ( l76l) . provides us with an analytic expression for r{j]) that we 
can invert to find rj (and hence f) as a function of r. 

Repeated differentiations of Eq. ( IT6l) with respect to r result in 

(dcj)-/dr])(dr]/dT) = iar, (A5a) 



(d(p-/dri)(d 2 ri/dT 2 ) + d 2 (j)-/dr l 2 )(dr l /dT) 2 = ia, (A5b) 

and so on, which when evaluated at S (where r\ = r = 0), yield dr]/dr\s, d 2 T]/dr 2 \s, etc., in 
terms of the known derivatives d<j)_/dr]\s, d 2 4>^/dr] 2 \s, etc., that constitute the coefficients 
in Eq. ( 1A4I) . Using these derivatives of rj at S, we can therefore write down the Taylor 
expansion of rj in powers of r: 

r) = If + f 2 + |f 3 + 8f 4 + 2§i f 5 + • • • , (A6a) 

where 

77 = lR P sin 3 ^ P sec 2 6 P r], (A6b) 
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and 

f = | exp(-i7r/4)i?2 s i n 5 9 P sec 4 6 P r. (A6c) 

The dependence of f on r now follows from Eqs. ( 1A2I) and ( 1A6I) . 

According to Eqs. ([63]), (El]), (174"]) . flAll and (ES), the explicit form of the function F| Ca 
is given by 



^| ?= _ (us2+iT2)1/2 ^#p 5 esc 5 # P u,(£ 5 2 + ir 2 ) 1 / 2 



x exp[-i(/i0 c - 7r/2)]9(r/- 3 )/9r 
for Rp ^> 1. Insertion of this expression in Eq. (ITU]) yields 

/ d£F(£, z) exp(i< 2 ) = §i?p 8 cot 6 fl P esc 5 # P exp[-i(/i0 s + tt/4)] 
</c 2 

x Uj(3Zi - T 2 ) 



(A7) 

c 2 



(A8) 



lev 



in which 



and 



/•oo 

/ drr" 4 exp(-ar 2 )^(r), (A9) 
Jo 



/■oo 

X 2 = / drr" 3 exp(-«r 2 )(d^/dr), (A10) 
Jo 

with ^ = (fj/f)^ 3 . Here, we have used the fact that /z0c — onis 2 = fufis, where 0s stands for 
the value of 0_ at S. 

That the integrals X\ and Z 2 have turned out to diverge is a consequence of our having 
interchanged the orders of integration and differentiation in Eq. ( |37j) [see also Eq. ( TlBl) ]. The 
standard technique for regularizing such divergent integrals is to treat them as generalized 
functions whose physically significant values (i.e., the values that we would have found 
had we not interchanged the orders of integration and differentiation) are given by their 
Hadamard finite parts 19]. 



To apply the technique to X\, one begins by appealing to Taylor's Theorem to represent 
the continuously differentiate function ip as 

^( r ) = ^(0) + V/(0)r + ±V"(0)r 2 + ^"'(0)t 3 + ^f"'( K r)r 4 , (All) 

where k is a number lying between and 1. One then inserts Eq. (1A11I) in Eq. (1A9I) to 
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rewrite X\ as 

/OO f'OG 
drr- 4 exp(-ar 2 ) +^'(0) J dr r~ 3 exp(-ar 2 ) 

+ ^"(0)X°°drr- 2 exp(-ar 2 ) + ^'"(0) J e °° dr r" 1 exp(-ar 2 ) ( A12 ) 
+ jr J e °° drij""(KT) exp(-ar 2 



The first four integrals inside the square brackets in this expression can be easily evaluated 
as functions of (a, e); e.g., 

/OO 
drr- 4 exp(-ar 2 ) = § e - 3 (l - 2ae 2 ) exp(-ae 2 ) + §7^3/2^^1/2^ (A13) 

in which the error function erfc(a 1 / 2 e) approaches unity in the limit e — > 0. 
The remaining fifth integral on the right-hand side of Eq. I I A 121) equals 



(A14) 



dr^ ,w (fi;r)exp(-ar 2 ) =4! J dr r" 4 [^(r) - ^(0) - i/j'(0)t - ^"(0)r 2 

- ^//"(0)r 3 ] exp(-ar 2 ) 



by virtue of Eq. (1A11I) . For a 1 (i.e., Rp ^> 1) and e = 0, the leading term in the 
asymptotic value of the right-hand integral in Eq. (1A14|) is given by 

POO 

drij""(KT)exp(-ar 2 ) ~ V""(0) / drexp(-ar 2 ) = Uit/a) 1 ' 2 ^""^). (A15) 



Here, we have applied l'Hopital's rule to remove the indeterminacy in the value of the kernel 
of the right-hand integral in Eq. (fill) at r = 0. 

Hadamard's finite part of the limiting version of each of the integrals that appear in- 
side the square brackets in Eq. ( 1A12I) is obtained by simply discarding those terms in its 
representation as a function of (a, e) that diverge when e tends to zero; e.g., 

JF| J™ dr r~ 4 exp(-ar 2 ) } = f Tr 1 / 2 ^ 3 / 2 (A16) 

according to Eq. ( 1A13I) . Thus, Eq. ( 1A15I) and the finite parts of the divergent integrals on 
the right-hand side of Eq. (IA12I) jointly yield 

T{Z X } ~f Tr 1 / 2 ^)^ 3 / 2 + !V'(0)(lna + 7 )a - |tt 1 / 2 ^"(0)« 1 / 2 
- i^"'(0)(lna + 7) + ^vr 1 /2^'"(0) a - 1 /2 ) a > i, 
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where 7 = 0.57721 is Euler's constant. The same procedure, when applied to the integral 
defined in Eq. (1A10I) . results in 

JF{J 2 } ~±V'(0)(lna + i)a - it 1 / 2 ^" {<S)a 1 ' 2 - ±^"'(0)(lna + 7) 

(A18) 

+ ^/^""{^a- 1 / 2 , a > 1. 
The required derivatives of i/j at r = can be read off the following expansion of (fj/r)^ 3 : 

i> = 1 - 3f - §f 2 - 4f 3 - iff 4 + • ■ ■ , (A19) 



which follows from Eq. (1A6al) [see also Eq. (1A6cj) . 



Evaluating the right-hand sides of Eqs. (IA17I) and flA18j) with the aid of Eqs. (lA6cl) and 



(1A19|) . and inserting the resulting expressions in Eq. (1A8I) . we finally arrive at 

^{3X1 - J 2 } = -(lOS/e^sin^eplsec^pl 15 ^ 572 , R P > 1, (A20) 
and hence, at Eq. (1771) . 
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